Quantum Dot Potentials: Symanzik Scaling, Resurgent Expansions and Quantum Dynamics 



Andrey Surzhykov, 1 Michael Lubasch, 1 Jean Zinn-Justin, 2 and Ulrich D. Jentschura* • 1 

1 Max-Planck— Institut fiir Kernphysik, Saupferchechveg 1, 69117 Heidelberg, Germany 
2 DAPNIA, Commissariat a VEnergie Atomique, Centre de Saclay, 91191 Gif-Sur-Yvette, France 

(Dated: February 6, 2008) 

This article is concerned with a special class of the "double-well-like" potentials that occur naturally in the 
analysis of finite quantum systems. Special attention is paid, in particular, to the so-called Fokker-Planck po- 
tential, which has a particular property: the perturbation series for the ground-state energy vanishes to all orders 
in the coupling parameter, but the actual ground-state energy is positive and dominated by instanton configura- 
tions of the form exp(— a/g), where a is the instanton action. The instanton effects are most naturally taken into 
account within the modified Bohr-Sommerfeld quantization conditions whose expansion leads to the general- 
ized perturbative expansions (so-called resurgent expansions) for the energy eigenvalues of the Fokker-Planck 
potential. Until now, these resurgent expansions have been mainly applied for small values of coupling param- 
eter g, while much less attention has been paid to the strong-coupling regime. In this contribution, we compare 
the energy values, obtained by directly resumming generalized Bohr-Sommerfeld quantization conditions, to 
the strong-coupling expansion, for which we determine the first few expansion coefficients in powers of g~ 2 ^ 3 . 
Detailed calculations are performed for a wide range of coupling parameters g and indicate a considerable over- 
lap between the regions of validity of the weak-coupling resurgent series and of the strong-coupling expansion. 
Apart from the analysis of the energy spectrum of the Fokker-Planck Hamiltonian, we also briefly discuss the 
computation of its eigenfunctions. These eigenfunctions may be utilized for the numerical integration of the 
(single-particle) time-dependent Schrodinger equation and, hence, for studying the dynamical evolution of the 
wavepackets in the double-well-like potentials. 

PACS numbers: ll.15.Bt, ll.10.Jj, 68.65.Hb, 73.21.La, 03.67.-a, 03.67.Lx, 85.25.Cp 



I. INTRODUCTION 

In this paper, we study the all-order summation of instanton 
contributions to the energy eigenvalues of anharmonic quan- 
tum mechanical oscillators which involve (almost) degenerate 
minima. The Euclidean path integral of quantum mechanical 
systems of this kind with one space- and one time-dimension 
is dominated by instanton configurations whose action re- 
mains finite in the limit of large Euclidean (imaginary) times. 
In order to find the energy eigenvalues, instanton configura- 
tions have to be taken into account. Modified quantization 
conditions have been conjectured for various classes of poten- 
tials (for a review see Refs.QlQ), and accurate numerical cal- 
culations have been verified against analytic expansions in the 
regime of small coupling^. Indeed, for small coupling, the 
energy eigenvalues are dominated by one-, two-, and three- 
instanton effects which correspond to trajectories of the clas- 
sical particle with a small number of oscillations between the 
(almost) degenerate minima of the potential. Here, we are 
concerned with the all-order summation of the instanton con- 
tributions, which is applicable to intermediate and strong cou- 
pling. Also, the transition from small to strong coupling, and 
strong-coupling expansions, will be discussed. Finally, we 
consider applications in the quantum dynamical simulation of 
finite systems. 

The purpose of this paper thus is threefold: first, to find 
generalized perturbative expansions (so-called "resurgent ex- 
pansions") for excited states of certain classes of notoriously 
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problematic 5 quantum mechanical potentials, second, to de- 
rive large-coupling asymptotics for these potentials and to in- 
vestigate overlap regions between small- and large-coupling 
asymptotics, and third, to outline applications of the consider- 
ations for the quantum dynamical simulation of a single parti- 
cle in a double-well-like potential. The first of these purposes 
is connected with mathematical physics, the second is tied to 
Symanzik scaling in the "poor man's" variant and therefore, 
to a basic implementation of the renormalization group, and 
the third one is rather application-oriented. 

Within the first and second aims of our investigation, 
we also investigate the fundamental question whether small- 
coupling perturbative expansions can be continued analyti- 
cally to the regime of large coupling, if one includes instan- 
tons into the formalism. Instantons can be considered either 
on the level of a resurgent expansion, augmented by an opti- 
mized (generalized) Borel-Pade resummation, or on the level 
of a generalized quantization condition, which entails an all- 
order resummation of the instanton expansion. 

The third application is mainly tied to the semiconductor 
"double quantum dot" structure a 6 ' 7,8,9 , which are formed from 
two quantum dots coupled by quantum mechanical tunnel- 
ing. Nowadays, these structures are generally accepted to be- 
long to one of the most hopeful candidates for the realization 
of quantum bits (qubits), because a single electron state in a 
double-well potential obviously can be localized in either of 
the two wells and, in that sense, represents a two-quantum- 
state system needed for quantum computing. Indeed, the the- 
oretical analysis of the structural and the dynamical properties 
of such (single-electron) double quantum dots can be traced 
back to double-well-like potentials. In this context, quan- 
tum dynamical calculations for a tunneling of a single par- 
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ticle between the two localized states nowadays attract special 
interest 10,11 but require a detailed knowledge on the eigenval- 
ues and eigenfunctions of the (double-well-potentials) poten- 
tials for a wide range of coupling parameters. 

Another question may as well be asked: The effective 
instanton-related expansion parameter, which reads Si (5) = 
y/2/ir c" 1 / 69 / ' g for the first two excited states of the Fokker- 
Planck potential (as discussed below), is nonperturbatively 
small for g — ► 0, but numerically not very small for some 
very moderate g. Specifically, Si (g) reaches its maximum 
Si(l/6) = 1.76115... already at a rather small coupling pa- 
rameter g = 1/6. So, one may ask how the "instanton expan- 
sion" in powers of Si(g) should be resummed, in addition to 
the perturbative expansions about each instanton. This latter 
step has never been accomplished, and we pursue its com- 
pletion via a direct resummation of generalized quantization 
conditions. 

This paper is organized as follows. In Sec. [Til basic defini- 
tions related to the Fokker-Planck and the double-well poten- 
tial are recalled. Calculations are described in SecHITI Specif- 
ically, we consider the resummation of the resurgent expan- 
sion in Sec. IIII Al the resummation of the quantization con- 
dition in Sec. IIII Bl large-coupling asymptotics in Sec. IIII CI 
and quantum dynamic simulations in Sec. HHP I Conclusions 
are drawn in Sec.lIVI 



H. BASIC FRAMEWORK AND NUMERICAL 
PROCEDURE 

A. Basic formulation 

In this manuscript, we discuss the determination of the 
eigenvalues of the one-dimension Fokker-Planck (FP) Hamil- 
tonian 
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where g is a positive coupling constant. For g = 0, Eq. 
represents the Hamiltonian of the quantum harmonic oscil- 
lator whose eigenvalues are given by the well known for- 
mula E^ K ^ = K, where K = 0,1,2,... is the "princi- 
pal" quantum number. For nonvanishing coupling, in contrast, 
no closed-form analytic expressions have been derived so far, 
and approximations have to be used (for a classification of the 
Fokker-Planck Hamiltonian in terms of a SUSY algebra, see 
App.lAl. The usefulness of the notation K instead of N will 
become clear in the following. If one considers the opera- 
tor V(g) = y/gq — ^/gq 3 + gq 4 /2 in 0} as a perturbation 
and formally applies Rayleigh-Schrodinger perturbative ex- 
pansion to the Kth harmonic oscillator state, then one finds 
the following result for the first terms, 



E, 



£ ) , P ert(9)=K-3K*g-[l7K* + ^K 



(2) 

All coefficients up to order g so are available for download 12 . 
This perturbation expansion^ fails to reproduce the spectrum 



of the Hamiltonian (Q]i even qualitatively. For instance, while 
the true ground-state energy E^" ' is manifestly nonvan- 
ishing and positive, the perturbation series (|2}, for K = 0, 
vanishes identically to all orders in the coupling g and is thus 
formally converging to a zero energy eigenvalue. A general- 
ization of perturbation theory is required, therefore, in order 
to correctly describe the physical properties of the Fokker- 
Planck Hamiltonian, including its energy spectrum. 

A complete description of the eigenvalues of the Hamilto- 
nian (fl} has been proposed recentlyi 2 ^ by using a generalized 
perturbation series involving instanton contributions. Since 
the concept of instantons in quantum mechanics has been pre- 
sented in a number of place a 1,2 ' 4 ' 13 , we may here restrict our- 
selves to a rather short account of the basic formulas. In 
the semi-classical framework, the eigenvalues of the Fokker- 
Planck Hamiltonian can be found by solving the generalized 
Bohr-Sommerfeld quantization conditional! 

1 
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In this expression, the functions Bpp(E,g) and App(E,g) 
determine the perturbative expansion and the perturbative ex- 
pansion about the instantons, correspondingly. 

The evaluation of these functions in terms of series in vari- 
ables E and g has been described in detail elsewher e) 1 ' 2 ' 13 
for rather general classes of potentials. In the particular 
case of the Fokker-Planck potential, for example, the func- 
tion Bpp(E, g) has the following expansion [see Eq. (14a) of 
Ref.H: 



Bpp (E, g) = E + 3E 2 g + 35E 



0(9 3 ).(4) 



The function Bpp (E, g) alone defines the perturbation expan- 
sion (f2]i which can be easily found by inverting the equation 
Bpp(E,g) = K. The instanton contributions to the eigen- 
values of the Fokker-Planck Hamiltonian are described by the 
function [see Eq. (14b) ofRef.H: 
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Extensive numerical checks of the generalized quantization 
condition (0 and the expansions <f4j — <j5j have recently been 
performed for the ground state of the Fokker-Planck potential 
in the weak coupling regime^. However, to the best of our 
knowledge, a numerical verification of these formulae (i) for 
excited states and (ii) for the large values of g is still missing. 
The numerical checks will be presented in Sec. IIII Bl 

While, of course, the present work is mainly devoted to the 
investigation of the energies and the corresponding wave func- 
tions of the Fokker-Planck potential, we will also briefly re- 
call the properties of the well-known double-well potential 
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FIG. 1: 
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(Color online) Eigenvalues E^p °' (left panel) and 
(right panel) of the Fokker-Planck Hamiltonian as a func- 



tion of the coupling parameter g. Results have been computed by the 
diagonalization of the Fokker-Planck Hamiltonian Q} in the basis of 
the harmonic oscillator wavefunctions. 



which is characterized by the Hamiltonian 



l -q 2 (1 - VOT) 2 • (6) 



Moreover, for the analysis of the energy spectra of the Hamil- 
tonians (Q]i and (|6]l it is very convenient to introduce the inter- 
polating potential: 

(7) 

which corresponds to the double-well potential if 77 = 0, 
whereas 77 = 1 gives the Fokker-Planck case. 



B. Numerical calculation of eigenenergies 

In order to numerically calculate energy eigenvalues of the 
Fokker-Planck and of the double-well potential, it is suffi- 
cient to consider matrix elements of these potentials in the 
basis of harmonic oscillator eigenfunctions, and to perform 
matrix diagonalization in a large basis spanned by harmonic 
oscillator eigenfunctions (typically, a number of > 5000 basis 
functions is sufficient for all calculations reported in the cur- 
rent article). One then observes the apparent convergence of 
the eigenenergies as the size of the basis is increased. This 
procedure is numerically stable provided one uses quadrupole 
precision (128-bit, 32 decimal figure) arithmetic. 

Fokker-Planck energies of the lowest three eigenstates [of 
the potential ([TJ] found by matrix diagonalization are dis- 
played in Fig. Q]as a function of the coupling g (we use the 
notation M = 0, 1, 2 in order to denote these three energy 
levels). As seen from this Figure, the states M = 1, 2 are de- 
generate in the limit g — * 0. A similar energy level splitting is 
well known for the symmetric double-well potential!^ and 
may be explained in terms of nonperturbative instanton con- 
tributions. In contrast to the double-well case, the Fokker- 
Planck potential contains a linear symmetry-breaking term 
[cf. Eq. Q)], and this term might be expected to lift any degen- 
eracy. However, excited states can still develop a degeneracy 
for g — > in view of the (only perturbatively broken) par- 
ity e = ±1 of the quantum eigenstates. The ground state of 
the Fokker-Planck potential, however, is located in one of the 



wells and does not develop any degeneracy due to parity (see 
also Fig.|9]below). 

It is interesting to investigate the adiabatic following of 
eigenvalues for the interpolating potential Q as a function 
of the parameter 77 for fixed g. This calculation (see Fig. |2]) 
reveals that the identification of the double- well energy eigen- 
valuesii^ w j tn quantum numbers (N, e) for the double-well 
with the quantum number M for the Fokker-Planck potential 
should proceed as follows: 



(N = 0,+)oM = 0, 
(N = 0,-)oM = 1, 
(N= 1,+)<^M = 2, 
(N= !,-)-&■ M = 3. 



(8a) 
(8b) 
(8c) 
(8d) 



The general relation is M = 2N + (1 — e)/2. However, 
the asymptotic behavior of the eigenenergies for g — > is 
different in the two cases: 



E. 



(N, S ) 
dw 

n(M) 



(g) -^N + 



1 

2 ' 

I(M + l)/2], 



(9a) 
(9b) 



where JxJ is the integral part of x, i.e., the largest integer m 
satisfying m < x. Equation d9bb implies that the perturbative 
contribution to the Fokker-Planck energy level with quantum 
number M is given by Eq. © with K = J(M + l)/2]. 

Apart from the degeneracy introduced by the instanton con- 
tributions, the eigenvalues E^>~ 1,2 \g) also have a quali- 
tatively different dependence on g when compared to the 
ground-state energy E^~ ^ (g). As seen from Fig.Q] while 
the energy E^~ ^ (g) increases monotonically as a function 

of the coupling constant g, the energies E^~ 1,2 \g) have 
minima at go w 0.07 and go ~ 0.025, respectively. 
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FIG. 2: (Color online) Adiabatic following of the lowest four eigen- 
values of the interpolating potential from 77 = (double-well) to 
77 = 1 (Fokker-Planck potential). The value of the coupling param- 
eter is held constant at g = 0.007. 
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III. RESUMMATIONS, ENERGY EIGENVALUES AND 
QUANTUM DYNAMICS 

A. Resummation of the resurgent expansion 

The generalized Bohr-Sommerfeld quantization condition 
(0 together with the expansions (Hji — (0 of the App and Bpp 
functions uniquely determines the eigenvalues of the Fokker- 
Planck Hamiltonian. For the ground state, the energy eigen- 
value can be found by systematic expansion of (O in powers 
of the two small parameters exp(— l/3g) and g, whereas for 
excited states, the parameters are exp(— l/6g) and g. Thus, 
for the Fokker-Planck potential, the particular form of the ex- 
pansion differs for the ground vs. excited states. 

Explicitly, the ground-state energy (M = 0) is given by the 
resurgent expansion^: 



E 



(9) = E 



n=l 



;-l/3fl 

2tt 
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k=0 



hi 



V f (0) a 1 

/ j J nkl 9 ' 



/=<) 



(10) 

where the index n denotes the order of the "instanton contribu- 
tion": n = 1 is a one-instanton, n = 2 is a two-instanton, etc. 
(as noted below, the one-instanton configuration involves a 
back-tunneling of the particle to the lower well for the ground 
state and thus has twice the action of the characteristic one- 
instanton effect for excited states). Another subtle point which 
should be recalled here is that the leading one-instanton term 
involves a summation over all possible n-instanton config- 
urations but neglects instanton interactions'. As seen from 
Eq. ( [Tol l, the evaluation of the ground state energy within such 
a (first-order) approximation, also referred as a "dilute instan- 
ton gas" approximation, requires the knowledge of the 
coefficients. Since these coefficients are available for down- 
load 12 , we only recall the six leading onesi 
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1 fW--- 
1 i J 101 — g ' 
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186624 ' 
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J 102 
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155 

3924815 
31104 ' 

163968231175 



(11) 



6718464 



By inserting these coefficients in Eq. (TTOb , we are able to 
perform now the numerical check of the validity of the one- 
instanton expansion for the ground state at small coupling. 
In Table U for example, the energy E^, (g) is displayed for 
coupling parameters in the range 0.005 < g < 0.03 and is 
compared to the "true" eigenvalues as obtained from the di- 
agonalization of the Fokker-Planck potential in the basis of 
harmonic oscillator wavefunctions. As seen from Table [H the 
ground state energy is dominated by the one-instanton effect 
for relatively small values of the coupling parameter g < 0.01. 
For stronger coupling, however, large discrepancies between 
the "true" energies and the results of resummation of Eq. ( TTOb 
at n = 1 are found indicating the importance of the higher- 
instanton terms which take into account the instanton interac- 
tions. The evaluation of the higher-order corrections (n > 2) 



TABLE I: Ground-state energy of the Fokker-Plank Hamiltonian. 
Results have been computed by the diagonalization of the Hamil- 
tonian in the basis of the harmonic oscillator wavefunctions ("true" 
energy) and by resummation of Eq. d 1 Ofc within the one-instanton ap- 
proximation (n = 1). The numerical uncertainty of the entries in the 
right column is estimated on the basis of the apparent convergence of 
results under an appropriate increase of the number of fjfy param- 
eters, which corresponds to the number of terms in the perturbative 
expansion about the leading instanton. Numerical discrepancies be- 
tween the left and right column are due to higher-order instanton 
contributions, as described in the text. We underline those decimal 
figures in the one-instanton results which are equal to the correspond- 
ing ones in the complete numerical solution. 



9 i?pp' °^ (diagonalization) n — 1 term of Eq. 1 1 Oi 

0.005 1.766 107 332 563 x 10~ 30 1.766 107 332 563 x 10" 30 

0.010 5.267 473 259 637 x 10" 16 5.267 473 259 637 x 10~ 16 

0.015 3.508 587 565 372 x lO^ 11 3.508 587 56 4 030 x IP' 11 

0.020 9.033 155 571 641 x 10 -09 9.033 15 4 730 920 x 10~ 09 

0.025 2.519 767 018 258 x 10~ 07 2.519 76 770 755(1) x 10~ 07 

0.030 2.313 302 179 961 x 10~ 06 2.313 251 574 075(2) x 10~ 06 



to the ground-state energy E F °p (g) is, however, a very diffi- 
cult task since it requires a double resummation of the resur- 
gent expansion, in powers of both g and exp(— l/3g). In the 
present work such a double summation based on sequentially 
adding higher-order instanton terms will not be performed. 
Still, for the sake of completeness, we here indicate the lead- 
ing two-instantor4 
Fokker-Planck ground state: 



AO) 

J 210 

f (0) 
J 211 



AO) 

J310 



AO) 

J311 



= 2. 



/ 2 ( o°o = 2 7 . 



2 ) J 2oi ~ 2 ' ° ■ 



= 87, .4o = 6 7 2 + 



(12) 



80 



= -f 7 - 6, = -15 7 2 - 12 7 - 17 + -7T 2 , 



where 7 = 0.577216 ... is Euler's constant. 

As seen from Eq. ( [Tol l, no splitting into levels with positive 
and negative parity arises for the ground state of the Fokker- 
Planck potential due to the linear symmetry-breaking term in 
Eq. (Q]i. This term modifies the potential in such a way that 
the leading, one-instanton (n = 1) shift of the ground state 
energy results from a back-tunneling (instanton-antiinstanton 
configuration) of the particle to the lower welli. For excited 
states, in contrast, the one-instanton configuration is a trajec- 
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tory which starts in one well and ends in the other, restoring 
the broken symmetry. Therefore, any excited state (M > 0) 
of the Fokker-Planck Hamiltonian can be characterized by its 
principal quantum number 



K = l(M + l)/2] 



and the parity 



2 2K — M 



(13) 



(14) 



In fact, this classification is very similar to the double-well po- 
tential © except, of course, the particular case of the ground 
state. It follows naturally that the resurgent expansion for 
the excited states of the Fokker-Planck potential is very close 
to the analogous expansion for the double-well potential and 
reads^: 



4p >0) (. 9 )-4p e) (3)-E^ 



19 



1=0 



71=1 



fc=0 



1=0 



where Eki are perturbative coefficients and z*k{9) is given 
by 



2^-1/2 



g K yJnK\ (K - 1)1 



-l/6g 



(16) 



The power of 5 can again be associated with the order of 
the instanton (K = 1: one-instanton, K = 2 means two- 
instanton, etc.). One should note that two intricacies are as- 
sociated to the precise meaning of the quantities that enter 
Eq. O: 

• In analogy to the double-well potential, the imaginary 
part which is generated by the resummation of the per- 
turbation series about the leading instanton (the "dis- 
continuity" of the distributional Borel sum in the termi- 
nology of Ref.ll4l) is compensated by an explicit imagi- 
nary part that stems from the two-instanton effect [from 
the factor ln(— 2/g)]. Related questions have been dis- 
cussed at length in Refs.QflU 

• In contrast to the ground-state energy ( fTOb , the leading 

contribution to the energies £&* >0) (s) far small cou- 
pling arises from the perturbation expansion (f2]i which 
is manifestly nonvanishing to all orders in g. However, 
since this perturbation expansion is independent of the 
parity e, the energy splitting of the levels with the same 
principal quantum number K is again dominated by the 
one-instanton contribution (n = 1). 

Similar to the ground state ( fTOb , we may compute such a con- 
tribution and, hence, a splitting of an arbitrary excited state K 



by making use of the fy^ coefficients, which for K > read 

(I = 0,1, 2, 3): 



f{K) 

J wo 



f(K) 
J 102 



AK) 

J 103 



=1, 

289 



hoi ~ 2 12 ' 

23 1139 45 
2 24 4 



(17a) 



695 
288 ' 



4913 
1885 



K" 



629 



K 



1637 



32 



A 



-K 6 



155825 
576 



-K — 



3835 
~2T 



-K — 



(17b) 



68885 

10368 ' 
(17c) 



Results for K < 28 are available for download^. In Table Ull 



1, 2 due to the one- 



for example, the splitting -E^p - 1 K ~ 
of the first two excited states M 
instanton effect (n = 1) is displayed as a function of the cou- 
pling parameter g. Again, a comparison of the results obtained 
by the resummation of Eq. (IT~5b and by the diagonalization of 
the Fokker-Planck Hamiltonian indicates the importance of 
the higher-instanton effects (n > 1) and, hence, the necessity 
of a double resummation of the resurgent expansion, in pow- 
ers of both g and exp(— l/6g). Instead of performing such 
a double summation explicitly, it is more convenient to enter 
directly into the quantization condition ((3), with resummed 
quantities as defined by the App(E, g) and Bpp(E, g) func- 
tions. We discuss this alternative approach in the next Section. 



B. Resummation of the quantization condition 

The resurgent expansions ( fTOb and (IT~5b for the energies of 
the ground and excited states of the Fokker-Planck Hamilto- 
nian follow as a direct consequence of the quantization condi- 
tion ([3j. As seen from our calculations summarized in Tables 
|I]and|lll these expansions are very useful for small coupling, 
but not of particular usefulness even for rather moderate val- 
ues of g, because of the necessity of their double resumma- 
tion. Here, we would like to investigate whether it is possible 
to resum the divergent series that gives rise to App(E, g) and 
Bpp(E, g) directly and look for solutions of the quantization 
condition (O without any intermediate recourse to the resur- 
gent expansion. In fact, this approach currently appears to 
be the only feasible way to evaluate the multi-instanton ex- 
pansion (in powers of n), because the quantization condition 
incorporates all instanton orders. 

In order to introduce such a "direct summation" ap- 
proach, we recall that the solution of the generalized Bohr- 
Sommerfeld quantization condition (0 for a particular cou- 
pling parameter g must provide the energy spectrum of the 
Fokker-Planck Hamiltonian. In other words, if one defines 
the left-hand side of the quantization condition (O as a func- 
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TABLE II: Energy difference between the excited states 
Epp~ 1,2 \g) of the Fokker-Planck Hamiltonian. Results have 
been computed by diagonalizing the Hamiltonian in the basis 
of the harmonic oscillator wavefunctions (left column) and by 
resummation of Eq. J 15b within the one-instanton approximation 
(right column). The numerical uncertainty is estimated on the basis 
of the apparent convergence of results under an appropriate increase 
of the number of parameters. As in Table[U we underline those 
decimal figures in the one-instanton results which are equal to the 
corresponding ones in the complete numerical solution. 



E ( F p~ 2) - -E*,"- 1 ' (diag.) n = 1 term of Eq. {15} 



0.005 9.848 553 978 903 x 10~ 13 

0.010 7.801 059 663 554 x 10~ 06 

0.015 1.213 924 539 483 x 10 -03 

0.020 1.289 613 568 640 x 10~ 02 

0.025 4.633 794 364 814 x 10~° 2 

0.030 9.699 341 140 782 x 10 -02 



9.848 553 978 903 x 10 -13 
7.801 059 65 9 99(1) x 10" OE 
1.213 9 1 452(2) x lO" 03 
1.288 765(1) x 10~ 02 
4^11 6(1) x 10~ 02 
9J>10 2(6) x 10~ 02 



tion of two variables E and g: 

Q{E l9 



T(-B FP (E,g))T(i- BME,g)) 



2\2B FP (E,g) e -A FP {E,g) 



fJ 



2tt 



(18) 



then the zeros of this function at fixed g determine the energy 
spectrum of (Q]): 

0(4^ (<?),<?) = 0, M = 0, 1, 2, (19) 

A numerical analysis of the function Q(E,g) can be used, 
therefore, in order to examine the validity and applicability of 
the generalized quantization condition given by Eqs. ^ — (0 
for the case of strong coupling. 

As seen from Eq. ( fT8l , any analysis of the function 
Q(E,g) can be traced back to the evaluation of the func- 
tions App(E,g) and Bpp(E,g) which constitute series in 
two variables, namely E and g [cf. Eqs. (|4]i — ©]. In or- 
der to compute these series, it is convenient to re-write the 
functions App(E, g) = App(E, gx)\ x= \ and Bpp(E,g) = 
Bpp(E, gx)\ x= i as (formal) power series in terms of a vari- 
able x, taken at x = 1 [cf. §8.5 of Ref. [Ell: 



N„ 



B FP (E,g) = ]T b$(E,g)x n 



n=0 



x=l 



A F p(E,g) 



39 



£ c$(E,g)x n 



x=l 



(20) 



(21) 




FIG. 3: (Color online) Energy dependence of (the real part of) the 
function Q(E,g). Calculations have been performed for fixed g — 
0.03. 



where the coefficients dpp (E, g) and bpp (E, g) are uniquely 
determined by Eqs. © and ©: byl(E,g) = E, 

b F p(E,g) = 3E 2 g, etc. In the computations, the power se- 
ries ( |20l and (f2TT > allow one to use a unified computer algebra 
routine for the Borel-like summations, which simply takes as 
input the variables a-pp (E, g) and bpp (E, g), as a function of 
g and returns the value of the resummed series at x — 1. In- 
deed, this routine can be universally used for different values 
of g and is therefore convenient for further numerical compu- 
tations which are discussed below. 

Making use of Eqs. < f2Qb and ( |2T1 i. we may now perform 
a simultaneous summation of the perturbation series as well 
as of the perturbation series about each of the instantons and 
find the functions Bpp(E,g) and App(E, g), correspond- 
ingly, where we use the same notation for a function and its 
Borel sum. There is a small subtlety because for positive g, 
the power series d20b and (I2TI 1 are nonalternating and diver- 
gent and, hence, special resummation techniques are required 
to calculate the Borel sums. In our present calculations, for 
example, we apply a generalized Borel-Pade metho d 16 ' 17 . The 
discussion of this method is beyond the scope of the present 
work, and the reader is r eferred to Refs. [l6llHl8l for a more 
detailed discussion. Because of their nonalternating prop- 
erty, the perturbation series defining the functions App (E, g) 
and Bpp (E, g) are Borel summable only in the distributional 
sensed. The evaluation of the Borel-Laplace integral thus 
requires an integration along a contour which is tilted with re- 
spect to the real axis (for details see Refs. [T7lll8lll9l and the 
contours C+i and C_i in Ref. [l6l) . The resummation of the 
divergent series (l20b and d2"TT i may be carried out along each 
of these contours, but it is important to characterize the pertur- 
bative and instanton contributions in the same way, i.e. to de- 
form the contours for Bpp(E, g) and App(E, g) either above 
or below the real axis, consistently. In the terminology of 
Ref. H one should exclusively use either "upper sums" or 
"lower sums," but mixed prescriptions are forbidden. From a 
historical perspective, it is interesting to remark that the possi- 
bility of deforming the Borel integration contour had already 
been anticipated in a remark near the end of Chap. 8 of the 
classic Ref. 

We are now in the position to analyze the properties of the 
function Q(E,g) and, hence, to extract the energy spectrum 
of the Fokker-Planck Hamiltonian. As mentioned above, to 
perform such an analysis for any particular g we have to (i) 



(») 



1 



TABLE III: Ground-state energy of the Fokker-Plank Hamiltonian. 
Results have been computed by the diagonalization of the Hamilto- 
nian in the basis of the harmonic oscillator wavefunctions (left col- 
umn) and by solving Eq. dl9t . as indicated in the right column. 



.'/ 


(diagonalization) 




ZeroofRe[Q(.E,g)] 


0.010 


5.267 473 259 637 x 


10" 


-16 


5.267 473 259 637 x 10~ 


0.030 


2.313 302 179 961 x 


10" 


-06 


2.313 302 17(2) x 10~ 06 


0.070 


1.267 755 797 982 x 


10" 


-03 


1.267 74(6) x 10~ 03 


0.100 


5.199 138 696 222 x 


10" 


-03 


5.199 3(2) x 10" 03 


0.170 


2.079 244 408 360 x 


10" 


-02 


2.078(1) x 10 -02 


0.300 


5.318 357 438 655 x 


10" 


-02 


5323(9) x 10" 02 



resum the (divergent) series for the functions App(E, g) and 
£?fp (E, g) and (ii) insert the resulting generalized Borel sums 
into Eq. ( fTSI l. We may then interpret the Q(E, g) as a func- 
tion of E (at fixed g) and (iii) numerically determine the zeros 
of this function which correspond to the energy values ELp 
of the Fokker-Plank Hamiltonian, according to Eq. dT9V In 
Fig. [3] for instance, we display the energy dependence of the 
real part of the function Q(E,g) taken at g = 0.03. In the 
energy range < E < 1, this function has three zeros which 
obviously correspond to the ground Epp and to the excited 

Epp 2 ^ states. The ground-state energy E F °p determined in 
such a way is presented in Table [III] and compared to ref- 
erence values obtained by the diagonalization of the Hamil- 
tonian matrix in the basis of harmonic oscillator wavefunc- 
tions. Moreover, apart from the particular case of g = 0.03, 
we also display the energy E F °p for other coupling parame- 
ters spanning the range from g = 0.01 to g = 0.3 (see also 
Fig- SI- This rather wide range of coupling parameters g con- 
sidered here allows us to investigate the behavior of the gen- 
eralized Bohr-Sommerfeld quantization condition in the tran- 
sition from weak to strong coupling. As seen from Table Hill 
the ground-state energy is well reproduced at g = 0.01 (up to 
14 decimal digits). Alternatively, a highly accurate value of 
the ground state energy (at g = 0.01) can be obtained from 
the the one-instanton contribution to the resurgent expansion 
([Tol l for g < 0.01, as indicated in Table Q] The accuracy 
of the one-instanton approximation is rapidly decreasing for 
higher g. For instance, at the moderate value of g = 0.03, 
the one-instanton term of the resurgent expansion ( fTOb repro- 
duces the ground-state energy only to four decimal digits (see 
the last row of Table [Q, while a total of eight digits can be 
obtained from solving Eq. dl9) . as indicated in the second 
row of Table [III] For even stronger coupling, one observes 
a much larger numerical uncertainty in the determination of 



g = 0.03 g = 0.1 g = 0.3 




2.0 2.2 2.4 2.6 4.5 5.0 5.5 4.5 5.0 5.5 6.0 
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FIG. 4: (Color online) Energy dependence of (the real part of) the 
function Q(E, g). As in Fig. [3] the results of linear regression anal- 
ysis of the function Q(E, g) are depicted by the solid line. Calcula- 
tions have been performed around the "true" energies of the ground 
state and for the different values of the coupling parameter: g = 0.03 
(left panel), g = 0.1 (middle panel) and g = 0.3 (right panel). The 
energy eigenvalues obtained using the displayed graphs determine 
the corresponding entries in the right column of Tablelllll 

the zeros of the function Q(E, g), because the convergence of 
the generalized and optimized Borel-Pade methods employed 
in the resummation of the App(E,g) and Bpp(E,g) func- 
tions is empirically observed to reach fundamental limits for 
larger values of the coupling, which cannot be overcome by 
the use of multiprecision arithmetic and might indicate a fun- 
damental limitation for the convergence of the transforms and 
are not due to numerical cancellations. It might be interest- 
ing to explore these limits also from a mathematical point of 
view. Specifically, we have determined the numerical uncer- 
tainty of the App(E, g) and Bpp(E, g) functions on the basis 
of the apparent convergence of the Borel-Pade approximants, 
inte grat ed in the complex plane and accelerated according to 
Ref.M7L using an optimal truncation of the order of the trans- 
forms. We found that as the order of the Borel-Pade transfor- 
mation was increased, the apparent convergence of the trans- 
forms stopped at around order 40 for g — 0.03 and higher. 
Despite these difficulties, the generalized Bohr-Sommerfeld 
quantization formula ( fl9l l determines the ground-state energy 
of the Fokker-Planck potential with an accuracy of about 0.01 
% up to g < 0.3 (cf. Table[m]). 

Until now we have discussed the computation of the 
ground-state energy Epp of the Fokker-Planck Hamiltonian. 
Of course, the function Q(E,g) may also help to determine 
the energies of excited states. In contrast to the ground state, 
however, the numerical analysis of the function Q(E,g) for 
excited states is more complicated due to bad convergence of 
the Borel sums for the App(E, g) and Bpp(E,g) functions 
in the energy range relevant for the excited states. As seen 
from Fig. [5] the convergence problems lead to relatively large 
numerical uncertainties for the numerical calculation of the 
function Q(E,g) already for a relatively mild coupling pa- 
rameter g = 0.07. We recall that this value of g corresponds 
to the minimum of the energy £pp /_1 ' ) (g) as a function of g 
and thus can be naturally identified as marking the transition 
from weak to strong coupling. As a result of the numerical 
uncertainties, the energy E^p^^g = 0.07) of the first ex- 
cited state may be reproduced only up to 2 decimal digits (see 
Table lFvl . For even larger values of parameter <?, the maximal 
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Energy [x 10 ] 



FIG. 5: (Color online) Energy dependence of (the real part of) the 
function Q(E,g). The results for the zeroes of Q(E,g) obtained 
by a quadratic regression analysis are depicted by solid lines. Cal- 
culations have been performed around the "true" energies of the first 
excited state with K = 1, e = +1 (M =1) and for three different values 
of the coupling parameter: g = 0.01 (left panel), g — 0.03 (middle 
panel) and g = 0.05 (right panel). The energy eigenvalues obtained 
using the displayed graphs determine the corresponding entries in the 
right column of TablellVI 



TABLE IV: The energy of the first excited state with M = 1 of 
the Fokker-Planck Hamiltonian. Results have been computed by the 
diagonalization of the Hamiltonian in the basis of the harmonic os- 
cillator wavefunctions (exact energy) and by solving Eq. < | 1 9b - 



9 (diagonalization) zero of the Re[Q(E, g)] 



0.010 9.677 074 461 352 x 10~ 01 9.677 074 461 352 x 10~ 01 

0.020 9.219 489 780 495 x 10~ 01 9.219 4 90(3) x lO" 01 

0.030 8.354 795 860 905 x 10~ 01 8.354 4(6) x 10~ 01 

0.070 6.828 548 309 058 x 10~ 01 6JS33(4) x 10~ 01 



accuracy of calculations, based on Eq. (T% , is only a single 
significant digit, even though the double resummation of the 
instanton expansion, and of the perturbative expansion about 
each instanton, is implicitly contained in the cited Equation. 

Supplementing the the specific Fokker-Planck energies of 
the M = 0, 1 states presented in TablesHHIand llVl we indicate 
in Fig. [6] the g-dependence of the energies 

EpP and £ F p 

as 

obtained from the analysis of the function Q(E,g) and com- 
pare to reference values obtained from the diagonalization of 
the Hamiltonian matrix. As indicated in Fig. [6] the results ob- 
tained by a direct resummation of the quantization condition 
give a good quantitative picture of the ground-state and the 
first excited-state eigenvalues in the ranges < g < 0.7 and 
< 9 < 0-3, respectively. The behavior of the numerical un- 
certainty as a function of g is explicitly shown in Fig. [71 where 
we plot the quantity 



A(ff) 



E 



(A/) 
FP. r. 



Afj) E- 



(M) 

FP, diag 



(g) 



E 



(M) 



FP, diai 



■(g) 



(22) 



with E 



(M) 



,(3) and 



-fp, rcsumvy; -fp, diag(ff) bein S the eigenvalues as 
obtained from the direct resummation of the quantization con- 
dition given in Eq. ( fT9l ) and the diagonalization of the Hamil- 
tonian (Q]), respectively (the latter values, which are numeri- 
cally more accurate, are taken as the reference values). While 
the accuracy of resummations for the ground-state energy re- 
mains satisfactory even for (relatively) strong coupling, the 
relative error for the first excited is numerically much more 
significant. 



C. Strong-coupling expansion 

As discussed in the previous Section, a numerical proce- 
dure based on the generalized Bohr-Sommerfeld formulae 
may provide relatively accurate estimates of the ground as 
well as the (first two) excited-state energies for the coupling 
parameters in the range < g < 0.3. The question arises 



0.250 9.508 936 793 119 x lO"" 1 9.4(2) x 10" 



whether g ss 0.3 can be considered as belonging to the strong 
coupling regime. Since the minima of the first two excited- 
state energies occur at g = 0.07 and g = 0.025, respectively, 
one might be tempted to answer the question affirmatively. 
However, one could devise a different criterion for the transi- 
tion to the strong coupling regime. For instance, one might 
define the strong-coupling region as a region of an appro- 
priately specified large-coupling asymptotic behavior of the 
eigenvalues Ep^ (g). 

The large-coupling asymptotics of the Fokker-Planck po- 
tential thus represents a natural next aim in the current inves- 
tigation. To this end, we apply a so-called Symanzik scaling 
q — > g~ x l e q in Eq. (Q]i and rewrite the Fokker-Planck poten- 
tial into another one with the same eigenvalues but a funda- 
mentally different structur e) 20 ! 21 : 



Hi 



FP 



tV3 



Ht, 



-1/3 



Q_-2/a 
2 y 



where the Hamiltonian does not depend on g: 



'I - 



(23) 



(24) 



We conclude that the Mth eigenvalue of the Fokker-Planck 
Hamiltonian for the g — * oo is determined in leading order by 
the Mth eigenvalue E^ I%> of the Hamiltonian (l24l : 



(25) 



Moreover, Eq. (l25T l also indicates that the classifications of 
the levels M = 0,1,2,... of the Fokker-Planck and the H h 
Hamiltonians are obviously identical in the strong-coupling 
regime. 

Based on Eq. d25l l. we now wish to compute the lead- 
ing asymptotics of the ground iJpp 1-0 '' (g) and excited-state 
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ground state M = 



first excited state M = ] 




0.2 0.3 0.4 0.5 
Coupling parameter g 




FIG. 6: (Color online) Eigenvalues E ] 



„(M=1) „(K=X,e= 
C FP — & FP 



-1) 



(M=0) 
FP 



(left panel) and 
(right panel) of the Fokker-Planck Hamil- 
tonian as a functions of the coupling parameter g. Results have been 
computed by the diagonalization of the Fokker-Planck Hamiltonian 
in the basis of the harmonic oscillator wavefunctions (solid lines) and 
by solving Eq. d!9t . where the latter correspond the the data points 
with the error bars. 



( AI— 1 21 

^fp \9) ener gi es of the Fokker-Planck Hamiltonian. 
This computation obviously requires information about the 
corresponding eigenvalues of the Hamiltonian H^. The en- 
ergies Ei M) have again been determined by a diagonalization 
of the Hamiltonian matrix within a basis of up to 1000 har- 
monic oscillator wavefunctions and than utilized in Eq. ( |25] >. 
As seen from Fig. [8j the leading asymptotics of the eigen- 
values Spp 7-0,1 ' 2 ^ (g) (dotted line), calculated in such a way, 
significantly overestimate the energies of the Fokker-Planck 
Hamiltonian for the region < g < 0.3. Higher-order cor- 
rections to the large-coupling asymptotics are therefore re- 
quired, in order to reproduce more accurately the asymptotics 
of the eigenvalues E F A p~°' 1,2 \g). We observe that we may 
apply standard Rayleigh-Schrodinger perturbation theory to 
Eq. d23l l and use the fact that the perturbative with respect to 

g V 2/3 /2, re- 



Hl, which is V(g) 



1/2 



'1/3 



mains Kato-bounded with respect to the unperturbed Hamil- 
tonian i?L for large g. Within such an approach, a strong- 
coupling perturbation expansion can be written for each en- 
ergy E^P (g): 



ff 1/3 E4 M V 

k=0 



■2fc/3 



(26) 



where L 



L 



(M) 
fc>0 



(M) 




E^ 1 ^ and the higher perturbation coefficients 



are calculated in the basis of the wavefunctions of the 

(M) 



unperturbed Hamiltonian d24l . The first six h\, ' coefficients 
for the ground M = as well as the first excited M = 1, 2 
states are given in Table [V] 

In Fig. [8] we implemented the first few expansion coef- 
ficients as listed in Table [V] in order to calculate strong- 
coupling asymptotics for the three lowest levels of the Fokker- 
Planck Hamiltonian (dashed lines). Obviously, the partial sum 
of the expansion d26*| i as defined by only six coefficients pro- 
vides much better agreement with numerically determined en- 
ergy eigenvalues than the leading asymptotics d25l l. For in- 
stance, the energy of the first excited state M = 1 is very 
well described by the (first six terms of the) strong-coupling 
expansion already at g = 0.07 (the agreement is better than 
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FIG. 7: (Color online) Relative error d22t of calculations of the eigen- 
values as a function of coupling constant g. Results are presented 
for the ground (M = 0) and the first excited (M =1) states of the 
Fokker-Planck potential. 



one percent). This point, which also corresponds to the min- 
imum of the function E^~ 1 \g), can thus naturally be iden- 
tified as the transition region between the regimes of weak 
and strong coupling. By using such a definition for the tran- 
sitory regime, we may finally conclude that the resummation 
of the Bohr-Sommerfeld quantization formulae ©-(O, us- 
ing the condition ( fT9] i, may provide reasonable estimates for 
the energy levels E^P (g) even in a limited subregion of the 
strong-coupling regime. 



D. Quantum dynamics 

In previous Sections of the current paper, we have presented 
a systematic study of the energy spectrum of the Fokker- 
Planck potential. In particular, two methods for computing 
the eigenvalues £7 P p (g) have been discussed in detail: (i) a 
"brute-force" method which is based on the diagonalization 
of the Hamiltonian matrix and (ii) a generalized perturbative 
approach which accounts for the instanton effects. While, of 
course, the precise computation of the energy levels is a very 
important task, the complete description of the properties of 
the particular Hamiltonian also requires an access to its wave- 
functions Q^p\q). Most naturally, these wavefunctions may 

be obtained together with the eigenvalues Ep^ ' (g) by matrix 
diagonalization. In our calculations, the basis of the standard 
harmonic oscillator wavefunction {4> n (q)}%Lo is used for the 



construction of a Hamiltonian matrix (<p n \Hpp\ 
Fokker-Planck eigenfunctions are thus given by: 



OO 
71=0 



<t>n(q) , 



The 



(27) 



AM) 



where the coefficients c n are found by the diagonalization 
procedure. In Fig. [9] we display, for example, the wavefunc- 
tions d27| i of the ground M = and the first excited M = 
1, 2 states as calculated for a coupling parameter g = 0.05. 



TABLE V: coefficients of the strong-coupling perturbation e 



pansion Q26) 
Hamiltonian. 



of the eigenvalues E F p °' 1 ' 2 " 1 of the Fokker-Plan 



M = 



M = 1 



M 



0.281067 805 
-0.132 985 313 

0.021367 333 
-0.000 876 935 
-0.000 060 335 
-0.000 001557 



1.854 587 292 
-0.209 853 650 
0.028174 214 
0.000 220 875 
0.000 031789 
0.000 000 739 



3.686 419 624 
-0.307985 031 

0.025 765 237 
-0.000 099 791 
-0.000 002135 
-0.000 000 385 



As is evident from Fig. [9] the symmetry-breaking term 
the Hamiltonian (Q]i leads to a ground-state wavefunctii 

^fp _ °' ) (<?) (dashed line) which is neither a symmetric nor a 
tisymmetric combination of the wavefunctions of the right ai 
left wells, but localized in the lower well. For the first excit 
states, in contrast, the symmetry is partially restored, and we 



may attribute the wavefunctions $p^f l ' (q) (dashed and dot- 
ted line) and $p^f -2 ' ) (q) (dotted line) to states with positive (e 
= +1) and negative (e = -1) parities [see also Eqs. (Qj) — (fl4b l. 



Until now, we have only discussed the evaluation of 
the eigenstates and eigenenergies of the Fokker-Planck 
Hamiltonian. In theoretical studies of double quantum-dot 
nanostructure a 10 ' 1 '- 22 i 23 and of the quantum tunneling phe- 
nomena in atomic physic o 24 ' 25 , a large number of problems 
arise, in which the time-propagation of some (specially pre- 
pared) wavepacket in double-well-like potentials has to be 
considered. The methods for such a time-propagation, such 
as the well-known Crank-Nicolson method, the split-operator 
technique, and approaches based on the Floquet formalism 
and many others, are discussed in detail in the literature (see, 
e.g., Refs. fl0ll22ll26ll27ll28b . As a supplement to our previous 
considerations, we will now consider the time-propagation 
of an initial wavepacket in a double-well-like potential, re- 
calling the adiabatic approach as one of the most simple and 
best-known techniques (see, e.g., the book|29| and references 
therein) for the integration of the (time-dependent) single par- 
ticle Schrodinger equation. Within such a technique, in which 
we can naturally make use of the results previously derived 
for the eigenstates and eigenenergies, the propagation of a 
wavepacket ^(q, t) in the (time-independent) Fokker-Planck 
potential ([]]) is given by: 



*(3,t) 



M=0 



6 (M) 



exp 



4^) 



(28) 



Here the coefficients b^ M ^ 



(M) 



^fp I = 0)/ determine 
the decomposition of the initial wavepacket (at t = 0) in the 
basis of the eigenfunctions (|27] |. 

Equation (l28l l provides an exact solution for the wavefunc- 
tion t) at an arbitrary moment of time only in the limit 
of an infinitely large basis of harmonic oscillator {(f> n (q)}^ =0 
and Fokker-Planck {$ F p f) (q)}f I=0 wavefunctions. For com- 
putational reasons, however, summations over basis functions 
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FIG. 8: (Color online) Exact values (solid line) for the ground state 
and the two lowest excited energy levels for the Fokker-Planck po- 
tential as a function of g, together with the leading asymptotics 
(dotted line) and the partial sum of the strong-coupling expansion 
(dashed line) as defined by the first six nonvanishing terms in powers 
of g~ 2fe / :! [see Eq. (|26}], which are listed in Table HITl Note that the 
leading strong-coupling asymptotics alone cannot satisfactorily de- 
scribe the true energy eigenvalues for moderate and small coupling. 
By contrast, the partial sum of the first six nonvanishing terms of 
the strong-coupling asymptotics yields numerical values which are 
indistinguishable from the true eigenvalues on the level of the line 
thickness used in the plots, down to rather small values of the cou- 
pling (dashed vs. solid lines). 



have to be restricted to a finite number. In our calculations, 
basis sets of 300-1000 wavefunctions have been applied de- 
pending on the parameters of the initial wavepacket and the 
coupling parameter g. The actual size of the basis has been 
chosen according to the numerical checks of the Ehrenfest 
theorem or the normalization of the wavepacket. 

Within the adiabatic approximation, which is valid for 
slowly varying potentials 29,3 ", we may divide the time evolu- 
tion of the potential into small intervals At and assume that 
for every fcth interval the Hamiltonian is time-independent 
and the propagation of the wavepacket is governed again by 
Eq. d28l l where, of course, eigenvalues E^^ and eigenfunc- 
tions $yp r ' ) (q) should be replaced with the eigenvalues Ef, 

and eigenfunctions (q) of the Hamiltonian Hk = H(tk)- 
We have applied this adiabatic time-propagation method, 
whose variations are well known from the literatur o 29 i 31 ' 32 ' 33 
and which is equivalent to an exponentiation of the instanta- 
neous Hamiltonian for each time interval At, to investigate 
the evolution of an (initially) Gaussian wavepacket in a time- 
dependent potential (0 which oscillates sinusoidally between 
the Fokker-Planck and the double-well cases. Since the ani- 
mated results of this simulation are available for download 12 , 
we just present a small series of snapshots in Fig.[l0] As seen 
from these pictures, the wavepacket, which is initially located 
in the right well, performs oscillations between the wells. 
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FIG. 9: (Color online) The wavefunctions $p P (q) of the ground 
M = and the first excited M = 1, 2 states of the Fokker-Planck 
Hamiltonian calculated at coupling parameter g = 0.05. The base 
lines for the plots of the wave functions correspond to their energies 



These oscillations are controlled by the temporal change of 
the potential 0. We have checked empirically that the adia- 
batic approximation employed here does not represent an ob- 
stacle for an accurate time evolution in even rapidly oscillat- 
ing potentials, because of the calculational efficiency of the 
other steps in our time propagation algorithm (notably, the di- 
agonalization including the determination of eigenvectors can 
be implemented in a computationally very favorable way on 
modern computers). It is thus possible to perform quantum 
dynamical simulations in potentials which oscillate between 
two limiting forms with two fundamentally different charac- 
teristic ground-state configurations, each of which is governed 
by instantons, though in a different way. A generalization of 
our approach to two-dimensional potentials appears to be fea- 
sible and is currently being studied. 

IV. CONCLUSIONS 

In this paper, we have investigated the Fokker-Planck 
[Eq. (0] and the double-well [Eq. (0] potential from the point 
of view of large-order perturbation theory (resurgent expan- 
sion and generalized quantization condition), in order to map 
out the regimes of validity of the instanton-related resurgent 
expansion for the lowest energy levels, and in order to ex- 
plore the possibility of reaching the strong-coupling regime 
via a direct resummation of the ^4pp and Bpp functions given 
in Eqs. (0 and (0, which enter the generalized quantization 
condition 0. The latter approach entails a complete double 
resummation of the resurgent expansions ( fTOb and ( fT5l ) both 
in powers of the instanton coupling exp(— 1/6</) and in pow- 
ers of the coupling g (perturbative expansion about each in- 
stanton). The quest has been to explore the applicability of 
resummed expansions for medium and large coupling param- 
eters, in the transitory regime to large coupling. 

It is quite natural to identify the transition region for the 
first excited state as defined by the minimum of the energy 
level E^~ 1 \g) as a function of g (see the right panel of 
Fig. [6j, which occurs near g w 0.07. As is evident from 
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FIG. 10: (Color online) Time propagation of the (initially) Gaussian 
wavepacket in the time-dependent potential Q with 77 = sin(t;/60) 
which oscillates between the double-well and Fokker-Planck poten- 
tials. The base line for the plot of the wave function corresponds to 
its instantaneous energy. 



Figs. [6] and it is possible to reach convergence for both the 
resummed secular equation $1% as well as convergence of the 
strong-coupling expansion d26| ) in a somewhat restricted over- 
lap region 0.04 < g < 0.3. The question whether it is possible 
to use the resummed instanton expansion for large coupling, 
cannot universally be answered affirmatively, although it is re- 
assuring to find at least a restricted region of overlap. In order 
to interpret the overlap, one should remember that for typi- 
cal Borel summable series as they originate in various con- 
texts in field theory, it is much easier to perform resummations 
at moderate and even large coupling parameters than for the 
instanton-related case considered here. An example is the 30- 
loop resummation of the anomalous dimension 7 function of 
six-dimensional 3 -theories and of Yukawa model theories 14 , 
which lead to excellent convergence for couplings as large as 
g = 10 and higher—. 

For the latter case, it is even possible to obtain, numeri- 
cally, the strong-coupling asymptotics on the basis of a re- 
summed weak-coupling perturbation theory (see, e.g., Ref.l35l 
for a remarkable realization of this idea in an extremely non- 
trivial context). The general notion is that any perturbation 
(a potential in the case of quantum mechanics and an inter- 
action Lagrangian in the case of field theory) determines the 
large-order behavior of perturbation series describing a spe- 
cific physical quantity. However, the potential or interaction 
Lagrangian also determine the large-coupling expansion for 
the physical quantity under investigation. This means that 
there is a connection between the large-coupling expansion 
and the large-order behavior of the perturbation series gener- 
ated in each theory, and this correspondence can be exploited 
in order to infer strong-coupling asymptotics even in cases 
where only a few perturbative terms are known^. Accord- 
ing to our numerical results, the corresponding calculation of 
strong-coupling asymptotics on the basis of only a few per- 
turbative terms would be much more difficult in cases where 
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instantons are present, even if like in our case, additional in- 
formation is present in the form of the generalized quantiza- 
tion condition (0). 

We deliberately refrain from speculating about further im- 
plications of this observation and continue with a summary 
of the application-oriented results gained in the current in- 
vestigation. Potentials of the double-well type are important 
for a number of application-oriented calculations, including 
Josephson junction qubit a 36-37 , inversion doubling in molecu- 
lar physic a 38 ' 39 , semiconductor double quantum dotsi2i2, as 
well as, in a wider context, Bose-Einstein condensates in 
multi-well trap a 40 ' 41 ' 42,43 . For the latter case, a number of 
theoretical works have been performed recently in order to 
explore the time propagation of (one particle) wavepackets in 
driven double-well potential a 10 ' 11 ' 23 . In Sec. Ill B I we discuss 
a numerical procedure for an accurate description of energy 
levels and of the corresponding wave functions, which can 
thus be used in order to construct basis sets for an accurate 
quantum dynamical time evolution of wave packets in both 
static as well as time-dependent potentials. This well-known 
adiabatic technique for the integration of the (single parti- 
cle) time-dependent Schrodinger equation is briefly recalled 
in Sec lIIIDl We illustrate this technique by a calculation for 
a potential which oscillates between the Fokker-Planck and 
the double-well cases, governed by a time-dependent interpo- 
lating parameter 77 = cos(uit) as given in Eq. (0. Similar 
calculations can be done for cases where the potential admits 
resonances as in the case of a cubic anharmonic oscillator. In 
this case, the method of complex scaling leads to a basis of 
states which can be used in order to start quantum dynamical 
simulations. Related work is currently in progress. 

V. ACKNOWLEDGMENTS 

U.D.J, acknowledges support from the Deutsche 
Forschungsgemeinschaft (Heisenberg program). M.L. is 
grateful to Max-Planck-Institute for Nuclear Physics for 
the stimulating atmosphere during a guest researcher ap- 
pointment, on the occasion of which the current work was 
completed. 

APPENDIX A: SUPERS YMMETRY AND THE 
FOKKER-PLANCK POTENTIAL 

This appendix is meant to provide a brief identification of 
the Fokker-Planck potential (Q~|) in terms of a supersymmetric 



(SUSY) algebra 4 ^. We define the operators 



with W(q) =q(l- ^fgq) and the SUSY Hamiltonian 



HsvsY - ( B B-B+ ) = { H H 2 ) • (A2) 
with 

Hi = K^) 2 + ^ 2 (i ~ ^ q)2 + ^ q -\' (A3a) 

H2 =~ \ (£) + \q 2 (1 - x/S9) 2 - Vgq+ \ ■ (A3b) 

Notice that W(q) finds a natural interpretation as a "superpo- 
tential" in the sense 

ofRefs.HS 

The Hamiltonians Hi and 
H2 are "superpartners." They are related to each other by a 
simple reflection and translation, q — ► l/s/g — q, and have 
the same spectra. In that sense, one may say that the Fokker- 
Planck potential is its own superpartner. The Fokker-Planck 
potential therefore constitutes a case of "broken supersymme- 
try" with zero Witten index [see, e.g., Eq. (2.88) of Ref. l45ll . 
The construction of the supersymmetric partner thus does not 
help in the analysis of the Fokker-Planck Hamiltonian. One 
is forced into the instanton-inspired analysis presented in the 
current study. 

As a last remark, we recall that the Fokker-Planck potential 



v FP (q) = (! - %/ot) 2 +Vot-2 = 2 - 

(A4) 

could be assumed to admit a zero eigenvalue. However, as 
is evident from the discussion following Eq. (7.28) of Ref.lH 
the corresponding eigenfunction is not normalizable, and thus 
cannot be interpreted as a physical state vector. Neither the 
Fokker-Planck potential nor its isospectral supersymmetric 
partner admit a zero eigenvalue. 
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